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A method of treating radiation problems in magneto-ionic (anisotropic) media is 
presented. A "wave matrix" is defined, the zeros of whose determinant are the propagation 
constants of the ordinary and the extraordinary plane waves. A derivation of the dyadic 
Green's function for the unbounded medium is given, which is also based on this matrix. 
A formula is arrived at, which gives the power radiated by any distribution of alternating 
current in terms of the wave matrix and the spatial Fourier transforms of the currents. 
The method is illustrated by a discussion of the power radiated by an elementary dipole. 



1. Introduction 

An ionized gas in a permanent magnetic field is an anisotropic dielectric medium. Two 
well known examples are the ionosphere and the plasma investigated in controlled thermo- 
nuclear fusion research. We shall call this medium a "magneto-ionic medium/' but other 
names, like "magneto-plasma," are also used in the literature. 

The dielectric properties of a magneto-ionic medium can be described by a dielectric 
permittivity tensor [1,2,3,4] 2 and the propagation of plane electromagnetic waves in such a 
medium has been subject to many invest igations [3, 5]. The purpose of this paper is to present 
a general treatment of electromagnetic radiation in magneto-ionic media. (The radiation 
properties of antennas and moving charged particles are, of course, modified by the anisotropic 
properties of the surrounding medium.) The method proposed avoids the introduction of 
vector potentials, Hertz vectors, anisotropic potentials, and the like. The computation of 
the fields excited by a known distribution of oscillating current will be reduced to elementary 
matrix operations and the evaluation of integrals. We are particularly interested in the 
(complex) power radiated by a current distribution and shall illustrate the method by discussing 
in detail the power radiated by an elementary dipole. 

We would also like to direct the reader's attention to other published methods |6, 7, 8, 9, 10, 
11, 12], most of which have been applied to special radiation problems, like Cerenkov Radiation, 
in magneto-ionic media. 

2. Background 

To simplify the analysis, we assume the medium to be homogeneous, of infinite extension, 
and nonmagnetic (Mrei=l)- All a-c quantities shall be described by their complex amplitude. 
If the factor exp (j<at) is dropped and MKS units are used, Maxwell's equations take the form 

VXE=-jco Mo H (la) 

VXH=jcoe eE+J, (lb) 



1 Contribution from Engineering Laboratory, 19 Parks Road, Oxford, England. 

2 Figures in brackets indicate the literature references at the end of this paper. 
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where J is the external current density — as produced by antennas or moving charged particles — 
which we assume to be known. The usual matrix formalism is used, and vectors are regarded 
as column matrices. The permittivity tensor is represented by the matrix e, whose components 
are in general complex numbers, to include conducting (i.e., lossy) media. For most applica- 
tions a matrix e of the simple form 



(2) 



can be used (see e.g., [2, 3]) . Here the 2-axis of the coordinate system is orientated in the direction 
of the applied permanent magnetic field. It can occur, however, that more complicated 
matrices have to be used, in order to describe the medium's physical behavior correctly [4]. 
We now rewrite eq (la) 

H=^vXE (3) 

a>Mo 
and eliminate the a-c magnetic field H from eq (lb) to obtain the wave equation 

^ w _ A l_^AE=-iwMcJ> (4) 

where 1 is the unit matrix, A=V 2 the Laplacian operator, and the dyade VV is standing for a 
matrix, whose elements are the differential operators &/dXidxjt. We have to solve this wave 
equation to obtain the amplitudes E (r) of the a-c electric field in every point of space produced 
by a given distribution of current. The solutions shall satisfy the condition that at great 
distances from the sources the fields represent divergent traveling waves. 

3. Plane Waves 

A set of simple solutions of the homogeneous wave equation (J=0) describe plane electro- 
magnetic waves, which we propose to discuss in this chapter. The planes of equal phase are 
specified by the wave normal 

n=(n h n 2 , n 3 ) 

whose components are the direction cosines 

Wi=sin a cos /3, 

7i 2 =sin a sin (3, 

n 3 = cos a. 

Let r=(x,y,z) = (xi, x 2 , £3) be the radius vector, drawn from the origin to any point in space, 
then 

n l x-\-n 2 y J rn z z=Yi'X= const (5) 

is the equation of a plane, and the amplitude of a plane wave with the wave normal n will vary 
as 

E(r)=E e- jk «' r . (6) 

With k being any complex number, it is convenient to define the vector 

k=kn. 
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We note that 

— e^ jk ' t =—jk i e' jk ' r ) (7) 

and, operating the matrices of the wave equation on the same exponential 

(w — Ai~ *) e- jk ' T =-i (k)e~ jk ' T , (8) 

we get the "wave matrix' ' 

i=kk-k 2 i+k 2 Q e., (9) 

Here is introduced the propagation constant of electromagnetic waves in vacuum 

The electric field of eq (6) has to satisfy the homogeneous wave equation. Using eq (8) and 
dropping the exponential, this condition can be written 

X(k)E =0. (10) 

As we look for nonvanishing fields, this condition can only be fulfilled if 

detX(k)=0. (11) 

Equation (11) determines the propagation constants of possible plane waves with given wave 
normal n. 

If e has the form as in eq (2), an elementary computation shows that 

det X (k)=ij(ei sin 2 a+e z cos 2 a) (k 2 -kj) (k 2 -kj z ), (12) 

with 



2 /y ^ 2 (^f — ej) sin 2 a + e 1 e 3 (l+ cos 2 a) ±-y/(e 2 1 — e 2 l — e 1 e 3 ) 2 sin 4 q + 4e 2 ! e 2 ; cos 2 a ( . 

k '<"l k «- 2(6, sin 2 « + 6 3 cos 2 a) "' {U) 

These relations show that in a magneto-ionic medium, as is well known, two types of 
plane waves are possible for a given wave normal. They are called the "ordinary" and the 
"extraordinary" wave, and their (complex) propagation constants are k r and k n respectively, 
whose values depend on the angle a between wave normal and permanent mag- 
netic field. The quantities k Ijn lk are known as the "refractive indices" of the corresponding 
waves. 

4. Dyadic Green's Function 

In this chapter we propose to present a solution of the wave eq (4) for any known current- 
distribution J(r) (confined to a finite region of space). Because of the linearity of Maxwell's 
equations there must be a linear relation between the components of a current element and the 
components of the electric field produced by the latter at a point r. We can, therefore write 
for the electric field E(r) produced by the entire distribution J(r') 

E(r)=fdr'#(r,r')J(r'), (14) 

+ 00 

where \dx f stands for \dx' dy' dz'. The matrix 6 is called "dyadic Green's function." 

— 00 

It has proved a powerful tool in treating other problems of electrodynamics and is well known 
for an unbounded isotropic medium [13, 14]. In the following we shall derive G for an un- 
bounded anistropic medium. An alternative derivation has been described by Bunkin [6]. 
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We use the following identities holding for Dirac's S-function 



r 



.+«. 



i-'^dku 



(15) 
(16) 



and the abbreviation 



8(r)=5(x)d(y)S(z). 
Because of eq (15) the electric field E of eq (14) satisfies the wave equation if (^satisfies 

(VV-A±-k£)G(r,r') = -ju> t i 18(r-r>). 



(17) 



We note that VV and A operate on the variables r only and not on r', and we have to assume 
that the interchange of these two operators with the integration dx' is permitted. 

A 

To find a suitable matrix G we multiply eq (8) by the inverse of the wave matrix from the 
right to get 



(VV— M—k 2 e)\- 1 e- jk ' r = — U- jk ' r . 



(18) 



We can do this for all real k\, k 2 , and k 3 , if we assume the medium to be at least slightly lossy. 
Then, the zeros of det X (k), kj, and kj If will have imaginary parts (i.e., plane waves are atten- 
uated). But as shown later the results are also valid for lossless media. 

We finally multiply eq (18) by e jk ' T ' and perform the integrations dk= I I dh\ dk 2 dk 3 to 



find that the matrix 



G^jdkX-^e-^- 



(19) 



satisfies eq (17). We have thus found the dyadic Green's function for the unbounded anisotropic 
medium. To this solution, of course, we can add any solution of the homogeneous wave equa- 
tion, if required by the behavior of the fields "at infinity. " It can be shown, however, that the 
result given in eq (18) satisfies the required condition that no incoming waves shall occur. 
The inverse of % can be computed by using Cramer's rule 



det X 



(20) 



where the matrix A is the "adjoint" of X. For a magneto-ionic medium with a dielectric tensor, 
as in eq (2), we find 



A(k)=k*n n-k*klL+lctE, 



(21) 



where £ is the adjoint of 



E= 



€l€ 3 , 


J€ 2 e 3 , 





— i e 2€ 3 , 


*l€3, 
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(22) 



and 

e i (n? +n$) + e 3 (nf+fll) ; ie 2 M +^i) + e 8 ?W e l n 1 n 3 +J€ 2 n 2 n 3 

A 

Z= — i€ 2 (7i?+7^) + €3^712; €!(w?+ri|) + €3(w|+r^); e l n 2 7h~je 2 n l n z (23) 

_ ^n{a z -^ 2 n 2 n z \ e 1 n 2 n z +je 2 n 1 n z ; ^(1+nf) 

The equivalence of 6r and Bunkin's Green's function can be shown by using eq (7) to replace 

each ki in A by the operation jd/bxt and interchanging the latter with the integration I dk. 

In this chapter the problem of finding the fields, produced by a know 11 distribution of cur- 
rent, has been reduced to the problem of evaluating integrals. The methods of integration will 
have to be adapted to the particular problem. For example, the method of steepest descents has 
been employed to find the dipole fields at great distances [6]. 

5. Power Radiated by a Distribution of Current 
The mean complex power radiated by the current distribution J(r) is given by 

P=-ifdrJ+(r).E(r), (24) 

where the row matrix J + is the Hermite conjugate of J. For some purposes it is convenient to 
rewrite this expression in terms of the wave matrix and the spatial^Fourier transform Jk of the 
current density 

J(r)=jdkJ k e~ jkr - (25) 

Using eqs (14), (19), and (25) we get 

P= ^J^J dk A # ft* dr dr > j+ x-i (k) J k „ 6 -;M*-k')-r'( k -k")] (26) 

With the help of the relations (15) and (16) these 15 integrals can be reduced to 3, and we 
obtain the relatively useful formula 

P=-4J7r 3 co M o jdk J+ X- 1 ^. (27) 

One can, of course, rewrite eq (24) in this form almost immediately, if Parseval's equation is 
used. 

An application of formula (27) is given in the next section. 

6. Power Radiated by an Elementary Dipole 

To simplify the problem, we assume here that the medium is lossless. Let an elementary 
electric dipole with moment p be placed at the origin of our coordinate system. The spatial 
distribution of current is then 

J(r)=jcop5(r). (28) 

For the calculation of many physical quantities, like the fields at great distances or the real 
power radiated, many current configurations J'(r), which are concentrated in an electrically 
small region, can be considered equivalent to an elementary dipole (see e.g. 115]). The 
equivalent moment is given by 

jcop=JrfrJ / (r). (29) 
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Certain quantities, however, like the fields in the immediate neighborhood or the reactive 
power, depend very strongly on the dimensions of this region. We are, therefore, mainly 
interested in the real power, radiated by the dipole. In an anisotropic medium the power is 
expected to be different for different orientations and polarizations of the dipole. 

As we want to apply formula (27), we need the Fourier transform of J(r), which is 



T 3<*V 



(30) 



Putting this into the formula, we see that the complex power can be written as a bilinear form 

(3D 



where the components of the matrix 






-stf**- 



(32) 



are measured in fim 2 (impedance units per unit area). This complex matrix z can be split 
into a Hermitian and an anti-Hermitian part 



z=r+jx, 



(33) 



that is to say, the matrices r and x are both Hermitian (r=r + , x=x + ). With the help of these 
two matrices the bilinear form of eq (31) splits up into two Hermitian forms (which are real 
numbers) thus separating real and reactive power. The real power P r is therefore given by 



r — j p+rp. 



(34) 



The reactive power produced by a region filled with current increases without limit as the 
region contracts. It has, therefore, no physical meaning for an elementary dipole. 

In the following we would like to sketch the steps of the computation of r for an e as in 
eq (2). That means, we have to pick the Hermitian part of I, given by eq (32), which is 

possible after performing two steps of the integration of dkX -1 . 

For this purpose we introduce polar coordinates in k-space, with the volume element 

dk=k 2 sin a dk da d/3. (34) 

We choose the intervals of integration from — oo to + oo for k, from to x/2 for a, and from 
to 2ir for j8, to cover all k-space. Equations (12) and (20) are used to substitute for 
X -1 , and eq (21) to rewrite 



k 2 A(k) 



=k 2 n n+(fcf+« r )n n-klt+~^- \jjz^ ^HP^fc ^ fe) } < 35) 



(k 2 -k 2 z )(k 2 -kh) -T W-T-W- -0- I k 2_ kh \J2- k 2 "W k 2_ ku 

As k T and k n are independent of f} 9 the integration with respect to this variable can be performed 
As a result of this, we introduce two new matrices, N and M, 

~sin 2 a 

1 C 2t 

nnJ/3= sin 2 a 

7T Jo 

2 cos 2 . 

*2€ 1 sin 2 a-f€ 3 (l+cos 2 a) ; 2^e 2 sin 2 a; 

— 2 i 7*e 2 sin 2 a; 2e 1 sin 2 a+e 3 (l-f cos 2 a); 

0; 0; 2e 1 (l + cos 2 o:)_ 
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a 1 f2 

N=- 

TT JO 



A 1 f 2tt a 

M=- Ld/3= 

7T Jo 



(36) 



(37) 



We have postulated a lossless medium. Therefore, e u e 2 , and e 3 have to be real numbers 

A A 

and the matrices M and N are Hermitian. kj and &fj are also real but & 7 and/or k n can be real 
or imaginary ("cut off" plane wave) because negative values of e x and/or e 3 can occur in mag- 
neto-ionic media. 

J + 00 
dk. From (35) we notice that the two 

last terms of the integrands have poles if ki and k n are real. As a consequence of this, the values 
of the corresponding integrals are not uniquely determined. And it is here that we have to 
remember the medium is regarded as at least slightly lossy (see sec. 4). The poles are then 
removed from the path of integration, and, we obtain with the help of Cauchy's formula and 
the residue concept 

dk 7tj 

'V 



r 



(38) 



k 2 -kf 

where the sign of the root of kj has been chosen such that Im & 7 < and Re ki^> 0, and 
similarly for k H . The result of the integration is now uniquely determined and can also be 
used for the lossless case. 

The first three terms of the integrand contribute to the anti-Hermitian part of z only, 
and can be dropped. Their integral does not exist (sec 4 remarks on reactive power). 

With w=cos a, the vacuum characteristic impedance Z = 
length \ =27r/Ar we obtain finally 



= Vmo/€o, and the vacuum wave 



. irZof 

r= — - .-=- I 



ilu 



F 



where 



2X8 Jo H[ tl +(e 3 -e 1 )u i }(k i I -kh) 



F=HeTm{(k 6 I -k 6 II )N-kl(k 3 I -kh)M+2(k I -k II )kfE} 



(39) 



The symbol "Herm" stands for "Hermitian part of." As A T , M, and E are Hermitian, this 
simply means that all terms with an imaginary k r or k, f have to be dropped. A*/ and k n vary 
with ^( = cos a) according to eq (13) and can be real for some regions of integration and imag- 
inary for others. The latter is the case for angles a between wave normal and permanent 
magnetic field where the corresponding plane wave is "cut-off" and cannot transmit power. 

A A A 

Because of the particular structure of N, M, and E the matrix /• has the form 



r= 



>1 


-j>'2 


0" 


J r 2 


r x 
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(40) 



from which the power P r , radiated by dipoles of any orientation and polarization, can be 
computed. 

To discuss two special cases of polarization let us consider: 

(a) A linearly polarized dipole moment p, as produced by an oscillating charged particle 
or an electrically small linear antenna. In this case the terms with r 2 cancel in the bilinear form 
and the dependence of P T on the orientation of the dipole can be computed from 



P=-^ 
^ T 2 P 







(41) 



which gives, if plotted in a polar diagram, an ellipsoid of revolution for l/-y/P T . For the special 
orientation perpendicular to the magnetic held the power radiated is 



551273— CO- 
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and for parallel orientation we find 



Pr \\—~2 r *\ 



(b) A circularly polarized moment, as produced by a charged particle on an electrically 
small circular orbit or two crossed linear antennas with a difference of tt/2 in phase. We 
propose to discuss two particular orientations only. The first is such that p z = 0, and p y = ±jp x 
(for left- or right-hand polarization respectively) where we obtain 

Pr=^\ Vx \\r l ±r 2 ) 

which is a characteristic result for magneto-ionic media. It shows that the two crossed 
antennas are "coupled" by the medium. The power they radiate is not the sum of the powers 
that each individual antenna would radiate, if the other one were not excited. 

The second orientation is one where p y =0 and p z = ±jp x . Here the total power 

CO 2 

p r=^ \PxWn-\-n) 

is equal to the sum of the powers that would be radiated by the individual antennas. 

The values of r h r 2 , and r 3 are determined by the corresponding integrals of eq (39), which 
have yet to be evaluated. For the special case of an isotropic medium (e 2 = 0, e 1 = e 3 ) 

2xZ / — 
72=0; n=r 3 =y ^- ^€ Q =r Q 

if e^es^O, and r 1 =r 2 =7 , 3 =0 if ei = e z <C0, which is the case in an isotropic plasma at frequen- 
cies below the plasma frequency. If e 2 = and e 1} e 3 >0, which occurs in magneto-ionic media 
with a very high gyrofrequency or in an uniaxial crystal, the integration yields 

ri/r =i(3+€ 3 /€ 1 ); 

n/r =l; r 2 =0. 

With the exception of these and a few other special cases, the integrals have to be evalu- 
ated numerically. Some plots of an evaluation of rjr and rjr by means of an electronic com- 
puter can be found in [16]. 

If the power radiated and the current distribution of an antenna is known, the radiation 
resistance (which is defined by them) can be easily computed. Take for example a short linear 
antenna of length I and with a constant distribution of current /, which is orientated perpen- 
dicular to the permanent magnetic field. It produces a moment with the components 
p x =Il/jco, Py=p z =0, and its radiation resistance is R=l 2 ri. 



The author thanks H. Motz for encouragement and valuable discussion. 

7. References 

[1] H. W. Nichols, J. C. Schelleng, B.S.T.J. 4, 215 (1925). 

[2] W. P. Allis; Motions of electrons and ions, Handb. Phys. 21„ 383 (1956). 

[3] J. A. Ratcliffe, The magneto-ionic theory and its application to the ionosphere (Cambridge Univ. Press, 

Cambridge, England, 1959) (with further references). 
[4] J. E. Drummond, Basic microwave properties of hot magneto-plasmas, Phys. Rev. 110, 293 (1958) (with 

further references). 
[5] M. Born, Optik, pp. 413-420 (Berlin, Germany, 1933). 

[6] F. V. Bunkin; On radiation in anisotropic media, J. Exptl. Theoret. Phys. (USSR) 32, 338 (1957). 
[7] A. A. Kolomenskii, Radiation from a plasma electron in uniform motion in a magnetic field, Doklady 

Akad. Nauk S.S.S.R. 106, 982 (1956) (with further references). 

522 



[8] A. G. Sitenko, A. A. Kolomenskii, Motion of a charged particle in an optically active anisotropic medium, 

J. Exptl. Theoret. Phys. (USSR) 30, 511 (1956). 
[9] K. A. Barsukov, On the Doppler effect in an anisotropic and gyrotropic medium, J. Exptl. Theoret. Phys. 

(USSR) 36, 1485 (1959). 
[10] G. A. Begiashvilc, E. V. Gedalin, Cerenkov radiation of a magnetic dipole in an anisotropic medium, 

J. Exptl. Theoret. Phys. (USSR) 36, 1939 (1959). 
[11] V. L. Ginzburg, V. Ya. Eidman, The radiation reaction in the motion of a charge in a medium, J. Exptl. 

Theoret. Phys. (USSR) 36, 1823 (1959). 
[12] A. D. Bresler, The far fields exited by a point source in a passive dissipationless anisotropic uniform wave- 
guide; Trans. I.R.E. MTT-7, 282 (1959). 
[13] F. E. Borgnis, Ch. H. Papas, Randwertprobleme der Mikrowellenphysik, p. 251-60 (Springer, Berlin, 

Germany, 1955) (with further references). 
[14] P. M. Morse, H. Feshbach, Methods of theoretical physics, Pt. II, Ch. 13 (McGraw-Hill Book Co., Inc. 

New York, N.Y., 1953). 
[15] J. A. Stratton, Electromagnetic theory, p. 432 (McGraw-Hill Book Co., Inc., New York, N.Y., 1941). 
[16] H. Kogelnik, The radiation resistance of an elementary dipole in anisotropic plasmas, Fourth Intern. 

Conf. on Ionization Phenomena in Gases; (Uppsala, Sweden (North-Holland Publishing Co., 

Amsterdam, Aug. 1960). 

(Paper 64D5-90) 



523 



